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We analyze the semiclassical evolution of Gaussian wavepackets in chaotic systems. We prove 
that after some short time a Gaussian wavepacket becomes a primitive WKB state. From then 
on, the state can be propagated using the standard TDWKB scheme. Complex trajectories are 
not necessary to account for the long-time propagation. The Wigner function of the evolving state 
develops the structure of a classical filament plus quantum oscillations, with phase and amplitude 
being determined by geometric properties of a classical manifold. 
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Introduction. The standard approach to semiclassi- 
cal evolution is time-dependent WKB theory (TDWKB) 
[11) S H) Bl ■ This theory provides a clear geometric de- 
scription of the dynamics: a time dependent quantum 
state is associated with an evolving Lagrangian mani- 
fold in classical phase space. [A phase-space manifold 
{p{q)jQ) is Lagrangian if p = \/S{q), for some generat- 
ing function 5 0].] For TDWKB to be applicable the 
initial state must itself be related to a Lagrangian man- 
ifold. This is the case, e.g., of eigenstates of position or 
momentum operators (related to planes) or highly 
excited eigenstates of bounded integrable Hamiltonians 
(related to tori) Q . 

Somewhat surprisingly TDWKB has never been used 
to propagate Gaussian wavepackets. It is probably the 
static point of view what blocked the use of standard 
TDWKB: If a Gaussian state is thought of as the ground 
state of some harmonic oscillator, then it certainly does 
not qualify as an initial WKB state. However, by tak- 
ing the dynamics into consideration, a new perspective 
arises. In chaotic systems, when one observes the evolu- 
tion of an initially Gaussian wavepacket through a phase 
space representation (like Husimi or Wigner) 0], it be- 
comes manifest that, after some time, the wavepacket 
acquires the form of a thin filament, very similar to the 
classical evolution of the initial density 0,[l,|^,[l3] (in the 
case of the Wigner function the filament is decorated by 
interference fringes). The smaller h (as compared with 
the relevant action scales), the stronger the localization 
of the wavepacket along some classical manifold. With 
this picture in mind it is natural to conjecture that the 
wavepacket evolves into a WKB state, its support being 
a real phase-space manifold d, [ill, [12] ■ The purpose of 
this paper is to prove this statement. 

We show that, after some (short) time, a Gaussian 
wavepacket becomes a primitive WKB state. From then 
on, the state can be propagated using the standard TD- 
WKB scheme. Complex trajectories [ij, [l3, [IB] are 
not necessary to describe the long-time propagation of 



wavepackets, but they may be used to describe the evo- 
lution during the initial stage. The present approach not 
only offers an intuitive geometric description of the evolv- 
ing state, but can be very accurate, as we demonstrate 
with a numerical example. 

We focus on the Wigner function for its ability to re- 
flect subjacent phase space structures [l^, and because 
it is the standard representation when dealing with de- 
coherence (arising from coupling to the environment [^) 
in semiclassical regimes. 

TDWKB approach. In order to eliminate some un- 
necessary complications in the general problem of semi- 
classical wavepacket propagation in chaotic systems we 
consider a one degree of freedom Hamiltonian H{p, q, t), 
the dependence on t being periodic with period r. Let 
us assume that the one-period map, Mr(j>{t),q{t)) — 
(p(t + r), (/(i + r)), has a hyperbolic fixed point at the ori- 
gin, and, without loss of generality, choose the q axis to 
coincide with the unstable subspace. At i = we launch 
a Gaussian wavepacket at the origin: 



V'o(9) = (2^fT2)-i/4exp(-gV4a2) 



(1) 



Note that we have preserved the essential ingredient of 
chaotic dynamics, i.e., the exponential stretching and 
folding of phase space manifolds. 

If the g-uncertainty (cr) of the initial wavepacket ([T]) is 
large enough Je.g., a classical length, a ^ fi^], then ([T]) is 
already a primitive WKB state, 



V'o(g) = AQ{q)exp[iSo{q)/h] , 



(2) 



given that both the amplitude il) and phase 5*0 (g) = 
vary slowly on the quantum scale 2]. The associated 
Lagrangian manifold is p = = dSo/dq. Accordingly, 
this state can be successfully evolved using the TDWKB 
scheme, as we showed in Ref. [T7| . 

Let us now analyze what happens when a small cir- 
cular wavepacket is launched at an unstable fixed point. 
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Numerical simulations show that the positive part of the 
Wigner function gets stretched along the unstable man- 
ifold. As this positive part bends, interference fringes 
appear. The picture is that of a positive (classical) thin 
filament decorated by an oscillatory pattern (see Fig. [T]). 
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FIG. 1: (Color online) Linear density plot of the exact Wigner 
function [nh W{p, q)] after six periods of evolution with the 
kicked harmonic oscillator (see the text for a description of 
the system). We also show displaced versions of (i) the initial 
state (a circular Gaussian at the origin, left-bottom), (ii) the 
Wigner function after two steps (top-right). 

The fact that the structures in Fig. [1] are very similar 
to those found in the Wigner functions related to WKB 
eigenstates [ll], leads us to enquire: does the Wigner 
function of Fig. [T] correspond to a WKB state, i.e., can 
we write 

^,(g)«^A('^)(g)exp[^5('')(g)/?.]? (3) 

Here z/ labels the different branches of the hypothetical 
Lagrangian manifold supporting the WKB state {p'-"^ — 
dS^-^^ /dq) 3. It is understood that the expression above 
must be valid during some time interval, during which 
amplitudes and phases evolve according to TDWKB the- 
ory [1, [1,13, Hi J i-e., amplitudes are convected by the clas- 
sical flow in g-space 

Mq') = A^{q)\dq'/dq\-^'\ (4) 

and the generating function, which solves the Hamilton- 
Jacobi equation, can be written as an integral over the 
classically evolved manifold Q: 

St{q')^St{q^)^ r pdq. (5) 
J go 

In the case of only one branch, i.e., before the clas- 
sical manifold folds, a numerical simulation is sufficient 
to answer the question posed above, as we show in the 



following. For this purpose we introduce the concrete 
model system that will serve as our test bench. This is 
the kicked harmonic oscillator (KHO) 

H{p,q,t) = — \--muj'^q'^+Kcos(kq)'y' 5(t-nT) . (6) 
2m 2 ^ — ' 

The parameters 'm^uj = k = l, K = 2 (in appropriate 
units [10]), and r = tt/{3llj) guarantee a large chaotic 
region around the hyperbolic fixed point located at the 
origin The unstable direction is close to the q axis. 
We considered an initial circular state {t = 0^) given by 
Eq. (P) with cr = 0.08 (corresponding to ^ = 0.0128). 
Some snapshots of its evolution (Wigner functions) are 
shown in Fig. [T] 
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FIG. 2: (Color online) Wavefunction amplitude \ipt{q)\ and 
phase derivative d4>t/dq (inset) at t/r = 0,1,2,3,4 (black, 
green, red, blue, yellow, respectively). Also shown is the un- 
stable manifold (inset, cyan). Due to parity symmetry we 
only plot the region q <Q. 

In Fig. [5] we show amplitude pt (q) and phase deriva- 
tive d(j)t/dq of the evolved state for short times [ip = 
pexp(z^/?i)]. Both quantities must be smooth on the 
quantum scale for the state to qualify as a primitive 
WKB state. We see that as time grows the amplitude 
gets smoother: at t = it was localized in a region of size 
0{'h^^'^), at i = 3 it has acquired the maximal (classical) 
width. The derivative of the phase, which is the can- 
didate to Lagrangian support of the WKB state, stays 
smooth, and close to the unstable manifold, until t « 3. 
For larger times, ^ > 4, the turning point is reached and 
a fold is born, giving rise to quantum oscillations in both 
p{q) and (t>{q). For smaller h a similar behavior is ob- 
served, the only difference being that it takes longer to 
reach the turning point (this time goes like log h) . Our 
claim is that, for small enough h, there is a time window 
(imin, ^max) where ijjt is to good accuracy (see below) a 
primitive WKB state, meaning that it can be propagated 
further on according to TDWKB. In our numerical ex- 
ample, we checked that the optimal time for starting the 
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TDWKB scheme is t = 2. At t = 3 the interference ef- 
fects of the turning point are ah'eady significant (even if 
not apparent in Fig. [2]). At t = 1 the wavepacket is still 
not wide enough. 

Incidentally, the inset of Fig. [2] shows that TDWKB 
performs badly if applied at t = 0. The initial mani- 
fold p = does not evolve classically: it jumps over the 
unstable direction, instead of staying on the same side. 

The next step is thus to evolve ijjt=2 using TDWKB 
and compare with the exact propagation. In order to 
unfold subjacent phase-space structures, the comparison 
will be made at the level of Wigner functions. 

Semiclassical Wigner function. The Wigner function 
of a general WKB state [Eq. ([3])] can be calculated analyt- 
ically, using the stationary phase method, in a way analo- 
gous to that followed by Berry and Balazs for the special 
case of semiclassical eigenstates of integrable Hamiltoni- 
ans 0. We start with the usual definition of Wigner 
function, specialized to our case, 

W{p, 9,0 = ^/ d^At {q + OAt [q - ^e'*^^^^^ , (7) 

where the phase is given by 0(f) — St{q + S,) — St{q — £,) — 
2p^. We assume that (p, q) is not far away from the mani- 
fold. In this case only one branch of St is relevant Q • The 
stationary phase condition reads pt{q+S,)+pt{q—0 = 2p, 
where pt{q) = dSt/dq. If the point x = {p,q) is on 
the concave side of the evolved manifold but not too far 
away, there are in general two solutions ±f*, defining 
two points on the manifold, and X-. These are the 
tips of a chord having x as midpoint (Fig. [3]). When 
the stationary points are not coalescing they give indi- 
vidual complex conjugate contributions to the integral. 
The corresponding Wigner function reads 



Wip,q,t) 



Ao{q+)Aaiq^ 



COS {A/h - 7r/4) 



(8) 



The phase structure was extracted literally from Ref . [1] , 
i.e., A is the (symplectic) area between the manifold and 
the chord. The amplitude is different from Berry and Bal- 
azs , as they considered a different initial density. Here 
^o('Z±) is the initial amplitude at the preimages of x± 
and v± = dx±/dqa denote tangent vectors at x±, their 
moduli representing local rates of expansion (dqo is the q 
component of a displacement along the initial manifold, 
at the preimage of x±). 

Equation ([5]) is a good semiclassical approximation, ex- 
cept in the vicinity caustic points, where \v-\. Av- \ =0, 
i.e., when tangent vectors at the tips of the chords are 
parallel (see Fig. [3]). At caustics stationary phase points 
coalesce and one must use transitional (or uniform) ap- 
proximations [ll] . If (p, q) is close enough to the classical 
manifold, we can obtain a crude transitional approxima- 
tion to Eq. ^ as follows: (i) Approximate the manifold 
by the quadratic curve p ^ p^ +P*{q — q*)+P* (q — 9*)^/2. 
This leads to a cubic phase 4>{£,)- (ii) Neglect variations 
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FIG. 3: (Color online) Classical manifold corresponding to 
the state depicted in Fig. [T] (red). This manifold is itself a 
caustic of the Wigner function. The "ghost" lines (blue) also 
belong to the caustic. The point x is the center of the chord 
with tips at X- and x+. A point x belongs to the caustic if 
the vectors tangent to the manifold at x- and x+ are parallel. 
The dashed vertical line at g = —2.0 indicates a special section 
for testing the semiclassical approximations (see Fig. |4]). 



of the amplitude, i.e., At{q) ~ At{q*). If we traverse the 
caustic along the line q = q*, and assuming > 0, then 
HEI: 



W{p, q,,t) 



{n'p:y/^' 



Ai 



2(p-p*) 



{n'p:v/^ 



(9) 



Ai standing for the Airy function. 

We are now in position for the final step in our pro- 
gram: the comparison of exact and semiclassical evolu- 
tions. Figure [1] displays the exact Wigner function after 
six steps of evolution. It is evident that its skeleton is 
the caustic of the manifold of Fig. [31 which was obtained 
by evolving during four steps the manifold (phase deriva- 
tive) associated to the exact ^t=2{q)- 

For a quantitative examination, in Fig. [?] we plotted 
the section g = — 2 of the exact Wigner function together 
with the semiclassical prediction, Eqs. ([H]) and Inside 
their respective domains of validity both approximations 
are excellent. This completes our argumentation. 

Concluding remarks. We reported for the first time the 
use of standard TDWKB for semiclassical propagation of 
wavepackets in chaotic systems. 

The key point is that chaotic dynamics provides the 
initial expansion that defines the appropriate Lagrangian 
manifold for starting the TDWKB scheme. In particu- 
lar we showed that localized states typically evolve into 
WKB states, and explained how to calculate amplitudes 
and phases explicitly. 
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FIG. 4: A section of the exact Wigner function W{p,q = 
— 2) of Fig. [1] (line) vs. WKB approximation. Full circles 
correspond to the standard stationary phase result [Eq. ((Sjl] 
and open circles to the transitional approximation [Eq. ([9])]. 



Our results provide a novel perspective for re- 
examining important previous work on long-time 
wavepacket propagation. Consider, for instance, the re- 
markable calculations of Tomsovic and Heller's, who used 
a multiple linearization scheme to obtain accurate au- 
tocorrelation functions for large times [7|. In the light 
of our findings, their scheme can now be understood as 
arising from the linearization of the WKB wavefunction 
'iptil) (which is globally valid) in the vicinity of a peri- 
odic point. The family of "homoclinic" intersections in 
Ref. [7] , essential for organizing the summation of recur- 



rences, corresponds, in the TWKB context, to the set 
of intersections between the Lagrangian manifold of the 
evolved state and the stable manifold of the fixed point. 
In this way, we expect the present paper will contribute 
to the ongoing debate about the timescale for the break- 
down of semiclassical propagation. Whether the break- 
time diverges like \oe h fill . Il2l | or, more plausibly, like 
some power of h 0, |22|| is a question still awaiting a 
definitive answer. 

One important feature of TDWKB is that it can be 
easily supplemented to accommodate decoherence effects. 
Let us briefly analyse the example of the Lindblad master 
equation corresponding to a chaotic Hamiltonian system 
coupled to a high temperature reservoir rio!,[2l (or, sim- 
ilarly, the non-selective weak continuous measurement of 
q and/or jD l23l . [2^). In the formalism of quantum tra- 
jectories [2j|, the quantum jumps associated to such an 
environment (or weak measurement scheme}, amount to 
random rigid translations in phase space [25|. The al- 
ternation of random translations with Hamiltonian evo- 
lution leads to a final density matrix represented by a 
weighted ensemble of pure WKB states, each one related 
to a particular history of random translations. The corre- 
sponding Wigner function is thus suggestively expressed 
as an average over filamentary Wigner functions like the 
one depicted in Fig. [TJ 
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